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ABSTRACT 

During galaxy-galaxy interactions, massive gas clouds can be injected into the intergalactic medium which 
in turn become gravitationally bound, collapse and form stars, star clusters or even dwarf galaxies. The objects 
resulting from this process are both "pristine", as they are forming their first generation of stars, and chemi- 
cally evolved because the metallicity inherited from their parent galaxies is high. Such characteristics make 
them particularly interesting laboratories to study star formation. After having investigated their star-forming 
properties, we use photospheric, nebular and dust modeling to analyze here their spectral energy distribution 
(SED) from the far-ultraviolet to the mid-infrared regime for a sample of 7 star-forming regions. Our analy- 
sis confirms that the intergalactic star forming regions in Stephan's Quintet, around Arp 105, and NGC 5291, 
appear devoid of stellar populations older than 10 9 years. We also find an excess of light in the near-infrared 
regime (from 2 fim to 4.5 //m) which cannot be attributed to stellar photospheric or nebular contributions. This 
excess is correlated with the star formation rate intensity suggesting that it is probably due to emission by very 
small grains fluctuating in temperature as well as the polycyclic aromatic hydrocarbons (PAH) line at 3.3 yum. 
Comparing the attenuation via the Balmer decrement to the mid-infrared emission allows us to check the relia- 
bility of the attenuation estimate. It suggests the presence of embedded star forming regions in NGC 5291 and 
NGC 7252. Overall the SED of star-forming regions in collision debris (and Tidal Dwarf Galaxies) resemble 
more that of dusty star-forming regions in galactic disks than to that of typical star-forming dwarf galaxies. 
Subject headings: 



1. INTRODUCTION 

Collision debris are the offsprings of close encounters be- 
tween two or more galaxies. Their characteristics such as 
their morphologies or contents can vary greatly depending on 
the parameters of the interaction (prograde/retrograde, rela- 
tive speed, mass ratio, galaxy type, impact parameter, etc.): 
they might appear as long, narrow tails, diffuse bridges, shells 
and rings (Schombert et al. 1990). The smallest objects in col- 
lision debris are the size of HII regions, many of which have 
been discovered recently between galaxies in nearby groups 
and clusters (Gerhard et al. 2002; Yoshida et al. 2002; Sakai 
et al. 2002; Weilbacher et al. 2003; Ryan- Weber et al. 2004; 
Mendes de Oliveira et al. 2004; Cortese et al. 2006; Werk 
et al. 2008, 2010). These HII regions, also known as Emission 
Line Dots for their compact appearance, have a relatively high 
metallicity, consistent with the idea that the gas in which they 
lie was somehow ejected from parent galaxies. Larger ob- 
jects are often found in the tails and filaments resulting from 
tidal interactions between (gas-rich) galaxies. In the 10 6 M Q 

Electronic address: boquien@astro.umass.edu 

1 University of Massachusetts, Department of Astronomy, LGRT-B 619E, 
Amherst, MA 01003, USA 

2 Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, 
IRFU/Service d'Astrophysique, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, 
France 

3 Laboratoire d'Astrophysique de Bordeaux, OASU, Universite Bordeaux 
1, UMR 5804. CNRS/INSU. B.P. 89, F-33270 Floirac, France 

4 Dept. de Fi'sica Teorica y del Cosmos, Universidad de Granada, Granada, 
Spain 

5 Department of Physics, University of Crete, GR-71003, Heraklion, 
Greece 

6 IESL/Foundation for Research & Technology-Hellas, GR-71110, Her- 
aklion, Greece 

7 Chercheur Associe, Observatoire de Paris, F-75014, Paris, France 

8 NASA Herschel Science Center, California Institute of Technology, 
Mail code 100-22, 770 South Wilson Ave., Pasadena, CA 91 125, USA 



to 10 7 M mass range, some super star clusters can be formed 
(Weilbacher et al. 2003; de Grijs et al. 2003; Lopez-Sanchez 
et al. 2004). The largest HII regions are typically found at the 
apparent tip of tidal tails, at more than 100 kpc from the parent 
galaxies, in maxima of the HI column density where we also 
detect large quantities of molecular gas (Braine et al. 2000, 
2001; Lisenfeld et al. 2002) and where enough collision de- 
bris are often available to build a Tidal Dwarf Galaxy (TDG, 
"an object which is a self-gravitating entity, formed out of the 
debris of a gravitational interaction", Due et al. 2000). All 
these regions have the exceptional characteristic of being in- 
tergalactic star-forming regions and they are embedded in HI 
clouds out of which the stars have formed recently. 

Across their full range of masses, those debris - especially 
the ones that are gas-rich and host star-forming regions - have 
generated an increasing interest in the astronomical commu- 
nity over the past two decades. Indeed, for some of them 
the lack of a substantial evolved stellar population originating 
from the parent galaxies - that is, they are practically a sin- 
gle stellar population - makes them attractive targets to study 
star formation outside galaxies (Boquien et al. 2009, for an as- 
sessment of their use as laboratories to study star formation) 
to address a variety of problems such as constraining the pa- 
rameters of star formation and they also give hope to constrain 
the Initial Mass Function. At the same time, their high metal- 
licity, inherited from their parent galaxies, makes the study 
of their interstellar medium - in particular the detection of 
molecular gas through CO emission and dust - easier than in 
metal deficient dwarf galaxies. 

Studies of collision debris have addressed their molecu- 
lar gas content (Braine et al. 2000, 2001; Lisenfeld et al. 
2002, 2004), their dust emission (Higdon et al. 2006; Smith 
et al. 2007; Boquien et al. 2009), star formation rates (Han- 
cock et al. 2007; Boquien et al. 2007; Smith et al. 2008; Bo- 
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quien et al. 2009), stellar populations (Weilbacher et al. 2000; 
Mendes de Oliveira et al. 2004; de Mello et al. 2008; Smith 
et al. 2008; Werk et al. 2008; Hancock et al. 2009; Sheen et al. 
2009), dynamical properties (Bournaud et al. 2003, 2004; Due 
et al. 2004; Mendes de Oliveira et al. 2006; Amram et al. 
2007) and also dark matter content (Bournaud et al. 2007; 
Milgrom 2007). 

The recent availability of images of interacting galaxies by 
GALEX, Spitzer and Akari allows us to construct the SED 
(spectral energy distribution) of collision debris, with data 
points ranging from the far-UV (ultraviolet) to the mid-IR (in- 
frared). Fitting these SED with theoretical models, one may 
in principle determine several key parameters, such as the star 
formation rate, the star formation history, the attenuation as 
well as the total stellar and dust mass (taking into account the 
PAH [polycyclic aromatic hydrocarbons], the silicates and the 
graphites). In this paper, we present the first consistent UV 
to IR modeling of a sample of star-forming collision debris. 
Seven regions selected from the sample of intergalactic star- 
forming regions by Boquien et al. (2009) are studied. 

A number of questions remain open regarding star form- 
ing regions in the intergalactic medium, that can be addressed 
through such SED modeling. Are the star formation rate es- 
timators usually used for spiral galaxies (Kennicutt 1998) re- 
liable in intergalactic star-forming regions? Are there really 
collision debris that are constituted of a single, young, stellar 
population, that is that lack a billion years old stellar compo- 
nent? If so, can such simple regions be used to constrain the 
IMF (initial mass function) in an easier way than in other envi- 
ronments? What is the origin of their strong near-IR emission 
(as measured in the Spitzer 3.6 fim and 4.5 fim bands) if not 
due to an evolved stellar population? 

The sample and data are presented in section 2; the model- 
ing of their spectral energy distribution is described in section 
3; we present the results in section 4; we discuss them in sec- 
tion 5 and finally we conclude in section 6. 

2. SAMPLE AND DATA 
2.1. Sample 

Here our aim is to gain insights on star formation in in- 
tergalactic star-forming regions forming in collision debris, 
by modeling their spectral energy distribution. To do so, we 
have assembled a sample of intergalactic star-forming regions 
in interacting systems. The sample is drawn from the cata- 
log of star forming regions in collision debris published by 
Boquien et al. (2009). The objects in this catalog are specif- 
ically selected to span over a large range of star formation 
rates (from 0.002 M Q yr _1 to 0.2 M yr _1 as measured in the 
ultraviolet) in order to probe any variation from that of inter- 
galactic HII regions to that of candidate TDG. All the regions 
selected in the present paper have been confirmed to be in 
the intergalactic medium thanks to observations (HI, CO, UV, 
optical, IR, spectroscopy) but also thanks to numerical simu- 
lations for a number of systems. Moreover, spectroscopic ob- 
servations show they they all have a relatively high metallicity 
(8.4 < 12 + log O/H < 8.8), confirming that they are made 
from pre-enriched gas, and thus most likely from collision de- 
bris. In addition, we have also restricted the sample to those 
having high enough signal to noise data points in the UV, op- 
tical and mid-infrared wavelength domains, necessary for the 
modeling of their SED. Finally we have also made sure that 
all regions are dominated by photoionization processes and 
not shocks as we focus on star formation. This was performed 



measuring line ratios in available spectra of HII regions and 
comparing with photoionization models. This was particu- 
larly crucial for Stephan's Quintet which contains some re- 
gions that are strongly affected by shocks. These selection 
criteria give us a sample of 7 objects that will be studied in 
this paper. 

All systems are presented in Figure 1 . The first two regions 
we have selected belong to Stephan's Quintet. The first one, 
named SQ-2 in Boquien et al. (2009) is located at the tip of a 
gaseous+stellar tidal tail in the south-eastern part of the sys- 
tem. The second one, named SQ-5 is located in the HI tail in 
the north-east of the system. It is a low mass compact region 
with blue optical colors, resembling a series of point sources 
with strong spectral line emission. Arp 105 is a bona fide in- 
teracting system between a spiral galaxy and an elliptical one. 
There are two prominent tidal tails. The one towards the north 
is a mix of stars and gas. The one towards the south is likely 
constituted only of gas which makes it particularly interest- 
ing, with Arp 105S a star forming region at its tip. Arp 245 is 
an interacting system with a large stellar and gaseous tidal tail 
towards the north. We have selected Arp 245N which is a star 
forming region at its extremity. NGC 5291 is a system con- 
stituted of a large collisional ring containing numerous star- 
forming regions. NGC 529 IN is the brightest star-forming 
region in the ring and the brightest intergalactic star-forming 
region known so far in the ultraviolet, which is why we have 
selected it. Well-known system and typical of the last stage 
of the Toomre (1977) merging sequence, NGC 7252 has two 
very extended tidal tails. We have selected NGC 7252NW, 
a massive star-forming region at the extremity of one of the 
tails. Finally we have selected VCC 2062 which is a proba- 
bly the most evolved intergalactic region known to date but 
still forming stars (Due et al. 2007). These regions enable 
us to probe a large range of SFR in the intergalactic medium 
allowing us to investigate changes of the dust SED and the 
contribution of the evolved stellar populations to the infrared 
emission as a function of SFR. 

2.2. Data: from the UV to the mid-IR 

• All ultraviolet GALEX and Ha CFHT/ESO data fluxes 
are published in Boquien et al. (2009). 

• Optical broad-band images of Arp 105, Arp 245, 
NGC 5291 and VCC 2062 are respectively published 
in Due & Mirabel (1994), Due et al. (2000), Due & 
Mirabel (1998) and Due et al. (2007). Optical B, V and 
R-band images of Stephan's Quintet and NGC 7252 
were obtained by our group using the 3.5m at Calar 
Alto observatory and the NTT at La Silla ESO obser- 
vatory. 

• Near-infrared images of NGC 7252 and Arp 245 were 
obtained with the Akari space telescope and are pub- 
lished here for the first time. The N2 filter centered at 
2.4 /mi was used. Further ground-based near-infrared, 
J, H and K images of Arp 105, Arp 245 and NGC 5291 
are published by Due & Mirabel (1994, 1998); Due 
et al. (2000). 

• Mid-infrared Spitzer/IRAC data fluxes in the 3.6, 4.5 
and 8 fim bands are published in Boquien et al. (2009). 
Spitzer observations performed with the Short-Low 
module of IRS between 5 and 15 fim are available 
for Arp 105S, Arp 245N (Boquien et al. 2009) and 



3 




Fig. 1. — True-color images of the selected systems. In addition, the neutral gas is overlayed in blue and star formation through FUV, Ha or the 8 fim, is 
overlayed in red. More specifically in Stephan's Quintet: VLA/HI (blue), Spitzer/IRAC 8 fim (red), HST+Calar Alto 3.5 m (BVR). Arp 245: VLA-C+B/HI 
(blue), GALEX/FUV (red), Dupont 2.54m (UVI). NGC 5291: VLA-B/HI (blue), GALEX/FUV (red), DSS+NTT (BVR). Arp 105: VLA-C+B/HI (blue), 
Spitzer/IRAC 8 yum (red), CFHT/12k (BVR). NGC 7252: VLA-D/HI (blue), CFHT/Ha (red), DSS+NTT (BVR). VCC 2062: VLA-C/HI (blue), GALEX/FUV 
(red), optical SDSS+NTT (BVR). North is up, East is left. 



NGC 529 IN (Higdon et al. 2006). A 24 image 
of Stephan's Quintet obtained with Spitzer/MIPS is 
presented in Cluver et al. (2010). 24 yum images of 
NGC 5291 and NGC 7252 were obtained as part of our 
Spitzer's GO program 50191 and are published here for 
the first time. The VCC 2062 24 //m image was ob- 
tained from the public Spitzer archives. 

For all these images, aperture photometry has been carried out 
following the procedures described in Boquien et al. (2009). 
The apertures correspond to the ones presented in the afore- 
mentioned article, slightly widened if necessary to enclose the 
24 //m flux. For convenience we list in Tab. 1 (Sec. 4) the 
fluxes that we use in the present paper. 

Optical spectroscopy is available for all selected star form- 
ing regions but one. The slits are more narrow than the aper- 
tures and focus on the brightest parts of each region. The 
metallicity and attenuation of the ionized gas are published 
in the aforementioned articles. We have summarized the ob- 
served parameters in Table 2, Sec. 4. 

3. SPECTRAL ENERGY DISTRIBUTION MODELING 

In this section we model the spectral energy distribu- 
tion of 7 intergalactic star forming regions using a photo- 
spheric+nebular model and a dust model. 

We divide the spectral energy distribution modeling into 
two steps. In the first step we use Pegase II to fit the con- 
tinuum and line emission from the stars and the gas using ob- 
servations ranging from the ultraviolet to the near-infrared. 
The output of this fit provides information on the excitation 
and intrinsic interstellar radiation field. In the second step we 
use those estimates to fit the mid-infrared observations using 
a dust model. In the process we also estimate the energy ab- 
sorbed and re-emitted by the dust at longer wavelengths. 

3.1. Ultraviolet to near-infrared spectral energy distribution 

modeling 



We combine spectroscopic observations of emission lines 
with the Pegase II spectral evolution code (Fioc & Rocca- 
Volmerange 1997) to model star forming regions in collision 
debris. 

We use the Salpeter (1955) IMF from 0. 1 M Q to 100 M as 
a baseline for this study unless specified otherwise. 

3.1.1. Star formation history 

It is known that gaseous disks in parent galaxies are often 
more extended than the stellar ones. As a result, depending on 
the parameters of the interaction, the collision debris may be 
devoid of any older pre-existing evolved stellar population, or 
on the contrary be dominated by old stars. We model this us- 
ing a two burst star formation history. The first burst models 
the disk of a spiral galaxy, that is the evolved stellar popula- 
tion that is stripped during an interaction: it assumes a 12x 10 9 
years old exponentially declining burst which has consumed 
90% of its original gas reservoir. On top of this evolved popu- 
lation we add a second burst which models the star formation 
episode taking place in the collision debris after 12x 10 9 years. 
As the photometric properties change very rapidly during the 
first few million years the burst is modeled by an exponen- 
tially decreasing SFR: SFR(t) - ax exp(-t/r), where a de- 
termines the intensity of the burst and r is the timescale. We 
follow this second burst from t = 10 6 years to t — 2 x 10 9 
years, t being here the time elapsed since the beginning of the 
star formation episode in collision debris. For this popula- 
tion, t varies between 10 6 years and 10 x 10 9 years to span 
a large range of possible star formation histories from nearly 
instantaneous to nearly constant. 

The quality of the determination of the relative contribution 
of both components will ascertain the accuracy of the deter- 
mination of the other parameters. The relative importance of 
the two populations is determined by the ratio r, the stellar 
mass produced by the second burst divided by the mass of the 
stars produced by the first burst (the evolved population), at 
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a given time t for a given timescale t. We span a range of 
ratios from log r = -2.75 to log r — 2 by steps of 0.25. The 
parameter r can be seen as a burst strength. 

A possible limitation is the imprecise modeling of TP-AGB 
(thermally pulsating asymptotic giant branch) stars by Pegase 
II . Maraston (2005) has shown that they are dominant for 
populations aged between 300 x 10 6 and 2 x 10 9 years. For 
our specific case however, the older population is much older 
and is dominated by red giant branch stars and the younger 
populations best fit is never older than 300 x 10 6 years old as 
we will show later. 

3.1.2. Nebular emission 

The nebular emission accounts for a significant fraction of 
the flux in the broad optical bands from young star forming 
regions (Anders & Fritze-v. Alvensleben 2003). To improve 
the standard nebular modeling of emission lines by Pegase II 

- which assumes a given set of ratios between emission lines 

- we use attenuation-corrected emission lines obtained from 
spectra. To calculate the modeled flux in each line, we scale 
the observed spectrum so that the observed Hfi flux matches 
the modeled one. This is calculated using the following re- 
lation: F (H/3) = 4.757 x 10 -20 x N Lyc x /, where F(HJ3) is 
the flux in W, Niy C is the number of Lyman continuum pho- 
tons and / the fraction of Lyman continuum photons ioniz- 
ing an hydrogen atom usually taken as 0.7 (Mezger 1978). 
When there is no spectrum available, most notably in the near- 
infrared we keep the default Pegase II lines ratios. The neb- 
ular continuum is also modeled using Pegase II prescriptions. 
The attenuation of the lines is performed separately in a sub- 
sequent step (Sect. 3.1.4). 

3.1.3. Metallicity 

The metallicity has a significant influence on the stellar 
spectrum. Indeed, a low metallicity star is hotter and emits 
more ionizing photons, affecting the strength of the emission 
lines and star formation tracing bands. The star forming re- 
gions selected in the sample have a relatively high metallic- 
ity as shown by spectroscopic observations from the litera- 
ture (Due & Mirabel 1994, 1998; Due et al. 2000, 2007, see 
Table 2). We have assumed that the metallicity of the stars 
and the gas is constant. While local enrichment may be par- 
ticularly important at low metallicity, the effect is weaker in 
collision debris because of their high metallicity. We have 
therefore set the metallicity to the value observed spectro- 
scopically. 

3.1.4. Extinction 

Even though Pegase II has provisions to handle the atten- 
uation, we have performed this separately. We have mod- 
eled the attenuation of the burst taking place in the collision 
debris using the starburst attenuation law (Calzetti 2001), a 
standard law for strongly star-forming regions. Another law 
such as the SMC (Small Magellanic Cloud) law has a much 
stronger slope in the ultraviolet for a given Ay. However, such 
an attenuation law is unlikely as star forming regions in col- 
lision debris have a metallicity significantly higher than the 
SMC one (see Table 2). For the evolved population which 
has been stripped from the parent galaxies, we rather use 
the LMC (Large Magellanic Cloud) attenuation law (Gordon 
et al. 2003) with a foreground screen geometry. This law ex- 
hibits the typical 220 nm bump found in a number of spiral 
galaxies which is suitable for the evolved stellar population 



originating from the parent galaxies. In addition, this popu- 
lation shows up mostly in the optical and in the near-infrared 
part of the SED where differences between dust attenuation 
laws are not crucial. Thus, the choice of a different law for 
this population would only have a very small impact. For each 
population, we span independently an attenuation range from 
Ay — to Ay — 2.9 by steps of 0.1 magnitude. 

The separation between the attenuation of both components 
allows a simpler interpretation of the results and avoids intro- 
ducing new parameters to characterize the distribution of the 
dust and of the different stellar populations. In addition as 
star-forming regions are generally more extinguished than re- 
gions dominated by an evolved stellar population and are also 
much more compact, the error should remain minimal. 

3.1.5. Minimization 

To fit the models to the observations, we determine which 
one has the smallest^ 2 : 

2 , \ V \yi-axmi(au---,a M )] 2 
X (au---,a M )= >, , (1) 

where aj are the values of the parameters of the model, M is 
the number of parameters of the model, yi are the observed 
fluxes with uncertainties <r ( (assumed random, independent 
and distributed according to a Gaussian law), N is the number 
of measurement points, m, are the fluxes of the model to test, 
and a the following normalization factor: 

yt X mi (flu /of 

a ~ N 2 9 — P ' > 

Yji =1 mi(a\,---,aM) /erf 

All models produced by Pegase II are normalized to 1 M of 
gas. This normalization factor allows us to directly scale the 
model to the observations. 

This minimization is only performed on stellar population 
sensitive bands from FUV to near-infrared. 

3.1.6. Summary of free parameters 

The models have 5 free parameters: 1) the age t of the burst 
of star formation in collision debris, ranging from t — 10 6 
years and t — 2 x 10 9 years, 2) the time constant r of this 
burst ranging from r = 10 6 years (nearly instantaneous) to 
r = 10 x 10 9 years (nearly continuous), 3) the value of the 
attenuation of the evolved population of galactic origin mod- 
eled with the LMC law (Gordon et al. 2003), 4) the value of 
the attenuation of the younger population modeled with the 
starburst attenuation law (Calzetti 2001), 5) the ratio r of the 
mass of gas transformed into stars between the young burst 
and the evolved stellar population ranging from log r = -2.75 
to log r — 2. The total number of models considered for each 
star forming region is up to 127 x 10 6 . When there is only 
one stellar population formed in the collision debris, there are 
only 3 free parameters. 

3.2. Mid-infrared spectral energy distribution modeling 

Coupled with dust models, mid- and far-infrared data place 
constraints on the dust (size distribution of the grains, ioniza- 
tion of the PAH, etc.). We carry out such a study in star form- 
ing regions in collision debris. To do so we compare with 
templates of galaxies of various types and with dust models. 
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3.2.1. Galaxy templates 

As a first step, we perform a simple comparison between 
the observed UV-to-mid-IR SED of the collision debris and 
SED of template galaxies. Such a comparison is legitimate 
as our regions have luminosities and sizes of typical dwarf 
galaxies (Boquien et al. 2009). Galliano et al. (2008a) mod- 
eled the SED of 35 nearby galaxies of different types, sizes, 
star formation activities, and ages with metallicities ranging 
from a; 1 /50 to ^ 2 Z . We scale each template to fit the UV- 
to-near-IR broadband observations and select the one which 
minimizes the^ 2 . This process allows us to compare the in- 
frared SED of our sources to those of galaxies with similar 
stellar spectra. 

3.2.2. Dust modeling 

The intergalactic star-forming regions presented in this pa- 
per have been observed only in the mid-infrared and short- 
ward. We have fluxes only up to 8 /mi or 24 /mi on a select 
number of systems. The absence of far-IR data makes the 
modeling of these objects extremely speculative. Indeed, the 
far-IR emission provides the most important and necessary 
constraints for the modeling of the dust SED of star form- 
ing region in collision debris. The far-IR emission by dust is 
dominated by grains in thermal equilibrium with the radiation 
field. Therefore, for a given set of grain optical properties, 
the far-IR SED gives direct access to the distribution of equi- 
librium temperatures of grains located mainly in dense PDR 
(photodissociation regions) and in the diffuse ISM (interstel- 
lar medium). Consequently, it is the most reliable tracer of 
the physical conditions that the bulk of the dust is exposed to. 
On the contrary, the mid-IR emission is dominated by grains 
stochastically heated. The extent of their temperature fluctua- 
tions is unrelated to the energy density of the general radiation 
field. 



Normalised Dust SEDs 
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Fig. 2. — Dust SED of grains exposed to ISRF (interstellar radiation field) 
of various intensities. The grain properties (abundances, size distributions 
and optical properties) are from the bare grain model with solar abundance 
constraints of Zubko et al. (2004). The ISRF has the Galactic shape (Mathis 
et al. 1983), and is scaled by the factor Go (Go = 1 for the diffuse Galactic 
ISM). Each SED is normalized to its bolometric luminosity. 

Fig. 2 illustrates this phenomenon. It shows the dust SED of 
the diffuse Galactic ISM exposed to ISRF of different intensi- 
ties. The dust is made of carbon and silicate grains of different 
sizes (3.5 A to 0.25 /mi) responsible for most of the contin- 
uum emission, and of PAH containing 10 to 1000 C atoms, 



mainly responsible for the mid-IR emission features. The in- 
tegrated intensity of the ISRF is quantified by the standard Go 
parameter, defined by: 

Jr»0.2 fim 
4nJ A dA 
_ 0.0912 im 

1.6 x 10~ 6 Wirr 2 

where J A is the monochromatic mean intensity of the ISRF 
which is determined by the fit of the photospheric+nebular 
model from FUV to near-infrared. Using this ISRF, we com- 
pute the temperatures of our dust composition following the 
method of Guhathakurta & Draine (1989). Increasing the in- 
tensity of the ISRF makes the equilibrium grain SED shift to 
shorter wavelengths (because of hotter temperatures), while 
the spectrum of PAH and very small grains remains almost 
identical, since they are excited by single photon events. From 
Fig. 2, it is obvious that in the absence of far-IR observations, 
the mean intensity of the radiation field is unconstrained. 

Unfortunately, we do not have any constraints in the far- 
infrared. Consequently we do not pretend to constrain the dis- 
tribution of equilibrium grain temperatures within the galaxy. 
In addition, without further information, the diagnostics that 
we can extract from the mid-IR are very limited. However, 
we can put some constraints on the dust properties by relying 
on the information provided by the attenuation estimate (sec- 
tion 3.1 .4). Indeed, by assuming that the global emission from 
the star forming region is isotropic, the energy absorbed in the 
UV and in the optical - calculated from the attenuation deter- 
mined by the photospheric+nebular fit - should be re-radiated 
in the IR. 

Another important point is that the ISRF illuminate the 
PAHs. Indeed, for a typical cloud, the radiation field becomes 
redder, when the optical depth increases. The shape of this ra- 
diation field will have an effect on the emissivity of the PAHs, 
which are stochastically heated, but will not have any effect 
on the equilibrium grains which are only sensitive to the to- 
tal absorbed power. In our case, we consider the unattenu- 
ated ISRF (optically thin case), since PAH emission is dom- 
inated by photodissociation regions and cirrus clouds, which 
are both optically thin. This assumption introduces an addi- 
tional uncertainty on the determination of the PAH-to-dust 
mass ratio that we will derive in the next paragraph. How- 
ever, this error should be small compared to the variations we 
probe here as we will see in Sec. 5.2.2. 

Thus, we can perform a simple SED fit, where the total 
power emitted by the dust is determined by the energy bal- 
ance, and where the PAH-to-dust mass ratio is constrained 
by the mid-IR fluxes. Indeed, in first approximation, oc 
Mpah x G{). The fact that PAH are stochastically heated guar- 
antees that they emit at the same wavelengths. Hence, only 
one point is needed to scale their emission. The total in- 
frared power follows a similar relation: L/« oc Md ust x Go- 
Therefore the 8 /mi-to-IR luminosity is proportional to the 
PAH-to-dust mass ratio: Lg/Lip oc MpAjj/Mdust- As shown 
in Fig. 2, the PAH-to-dust mass ratio, /pah, is approximately 
independent of Go. Such an estimate of /pah is strongly de- 
pendent on the accuracy of the attenuation estimate, since we 
use the attenuation to guess the level of the far-IR emission. 
It is therefore not extremely reliable, in itself. However, the 
comparison of this derived PAH-to-dust mass ratio to its ex- 
pected value from chemical evolution considerations is more 
interesting, since it provides an insight on the amount of em- 
bedded IR emission. In other words, the attenuation estimate 
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from the visible recombination lines is biased towards evolved 
star forming regions, where the molecular cloud covering fac- 
tor is low enough for the radiation to escape. On the con- 
trary, the attenuation towards very young embedded regions 
will not be traced by the Balmer decrement. In consequence, 
comparing the derived PAH-to-dust mass ratio to its expected 
value is a way to estimate the fraction of embedded infrared 
power. Another way to probe the presence of embedded star- 
forming regions would be to compare the estimate of the dust 
masses with the attenuation. Unfortunately, this would re- 
quire to make assumptions on the geometry of the dust and 
the stellar populations, which is beyond the scope of this pa- 
per. This also prevents us from using the 24 //m emission to 
trace embedded star-forming regions as it depends on many 
unconstrained physical parameters such as the compacity of 
star-forming regions and the composition of the grains. The 
PAHs are a better tracer as they do not exhibit systematic tem- 
perature variations and finally we do not have 24 pm obser- 
vastions for all systems. 

4. RESULTS 

Spectral energy distribution fits allow us, as we have shown 
earlier, to determine a number of physical parameters such 
as the age, the attenuation, the star formation history. We 
apply the modeling technique described earlier to our sample 
of star-forming collision debris. 

We present in Tab. 1 the list of fluxes, corrected for the fore- 
ground Galactic attenuation, we use in this paper. Some fluxes 
might very slightly differ from the ones published in Boquien 
et al. (2009) when the aperture has been modified. 

The spectroscopically derived oxygen abundance, EW(H/3) 
and A v are given in Table 2. 

The best photospheric and nebular fit parameters are given 
in Table 3. Information obtained from the dust models is pre- 
sented in Table 4. 

4.1. Stephan's Quintet 

Stephan's Quintet shows signs of interaction (Moles et al. 
1997) with prominent intergalactic star-forming regions asso- 
ciated and large reservoirs of atomic and molecular hydrogen 
(Yun et al. 1997; Lisenfeld et al. 2002; Williams et al. 2002; 
Appleton et al. 2006; Cluver et al. 2010) but also hot stripped 
gas (Sulentic et al. 1995). We have decided to study two re- 
gions of particular interest. The famous and prominent objects 
called SQ-A and SQ-B are not studied here: their high optical 
attenuation make any fit of their SED very uncertain. See Fig. 
5 in Boquien et al. (2009) for an image of the system and the 
selected regions: 

• The so-called "tip-tail" region (region 2 in Figure 5 in 
Boquien et al. (2009) and hereafter) is located at the 
tip of the stellar tidal tail that emanates from NGC 
7319. Its optical morphology and rather red color re- 
veal the likely presence of dust lanes, and of an impor- 
tant old stellar component likely tidally expelled from 
the nearby spiral. Its vicinity to SQ-B in the tidal tails 
suggest that they have close chemical properties. 

• Region 5, on the contrary, appears as a very blue 
and compact star forming complex located towards a 
gaseous tidal tail. Optical images do not show any ev- 
idence of an evolved stellar population which makes it 
particularly interesting. 



4.1.1. Stellar populations 

The near-infrared band observations that are available in the 
archives are unfortunately too shallow to detect the selected 
star forming regions. However, given the apparent importance 
on optical images of the old stellar population in region 2, we 
have initially set the ratio r, of the mass stars produced by the 
second burst divided by the mass of the stars produced by the 
first burst, to the rather high value of log r — - 1 . 

The free parameters are the star formation history and the 
attenuation of both populations. 

As shown in Figure 3, the best fit (able to reproduce the ob- 
servations with a probability of 98%) is obtained with an age 
of t = 85 x 10 6 years, no attenuation for the evolved popu- 
lation of galactic origin, an attenuation of 1 . 1 magnitudes for 
the one formed in the intergalactic medium and a time con- 
stant t — 37 X 10 6 years. The Ha flux is also reproduced 
relatively well. The range of possible attenuations however is 
large: from 0.7 to 1.4 magnitudes. 

Assuming an attenuation of 1.1 magnitude, all time con- 
stants are possible with a probability of at least 75% to re- 
produce the observations. For a short timescale, typically 
t < 10 x 10 5 years, the age is about t ^ 30 x 10 6 years. 
Conversely a nearly continuous star formation yields an age 
comprised between t - 100 x 10 6 years and t = 200 x 10 6 
years. 

We have to note that the presence of a larger evolved pop- 
ulation would redden the optical spectrum, which as a conse- 
quence would make the model-derived age of the population 
formed in the intergalactic medium younger. The ages we find 
should then rather be seen as upper limits. In case there were 
proportionally fewer evolved stars, the ages would be older. 

For the blue - almost pristine - region 5, we have set log r = 

2. 

The fit in Figure 4 shows that unsurprisingly the burst of star 
formation is young, typically less than t — 10 x 10 6 years with 
the attenuation comprised between 1.1 and 1.2 magnitudes, 
relatively close to what has been measured spectroscopically, 
for a probability to reproduce the observations over 50%. The 
best fit (reproducing the observations with a probability of 
55%) corresponds to an age of t = 2 X 10 6 years and an atten- 
uation of 1.2 magnitudes. A burst age of a few tens of million 
years cannot reproduce the observations in any case. The age 
being young, the influence of the timescale is very marginal. 

4.1.2. Dust 

In Figure 5 we present the fit of the SED of region 2 with 
its best template and dust model. 

The closest template to the observations is M51, a metal- 
rich relatively dusty spiral galaxy. Although essentially con- 
strained in the UV-to-optical regime, the match is very good 
up to the mid-IR. The lack of observations beyond 24 pm pre- 
vents us from constraining Go- 

The fit of the SED of region 5, with its best template and 
dust model is presented in Figure 6. 

The best template in the UV-optical regime corresponds to 
Mrk 930, a low metallicity starbursting galaxy. In the mid- 
IR however, the fluxes in the intergalactic region are much 
higher, indicating important dust emission (consistent with its 
high metallicity). For the dust model, we show only one fit 
for a value of Go that was free to vary. The reason is that 
the 24 pm flux provides a constraint that cannot be explained 
with a simple very small grain continuum, without altering 
their size distribution. This fit is purely speculative, and we 
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TABLE 1 
Fluxes 





FUV 


NUV 


B 


V 


R 


J 


II 


K 


N2 


3.6 


4.5 


8.0 


24 






fim 




















fj m 




Stephan's Quintet 2 
Slcplian's Quiniei 5 
Arp 105S 
Arp 245 N 
NGC 5291N 
NGC 7252NW 


28.2+1.9 
5.5 + 1.4 
72.0 + 3.7 
83.7 + 3.8 
365.3 + 35.6 
39.7+1.3 


36.9 + 2.3 
8.8+1.4 
134.8 + 8.1 
154.1 ±2.5 
445.2 + 31.5 
41.5 + 1.8 


125.6 + 34.3 

12.3 + 2.5 
251.3 + 23.4 
1397.6 ±59.0 
754.3 ± 64.9 
91.3 + 11.8 


149.1 ±31.0 

15.4 ± 1.2 
292.8 ± 37.5 
2224.2 + 38.2 
834.5 ±64.5 
91.8 ± 11.4 


194.2 + 30.6 
23.9+1.3 

376.7 + 82.1 
2827.6 + 52.5 
978.2 + 109.0 

107.0+12.3 


286.2 ± 82.1 
7686.6 ±2526.7 
477.7 ± 45.0 


265.6 ± 202.8 
7287.1 ±3084.3 
455.9 ± 105.0 


289.9 ± 111.8 
8212.1 ±4140.4 
548.0± 181.0 


5872.4 + 880.9 
132.4+11.9 


215.3 ±50.0 
39.1 + 12.8 
489.3 + 191.3 
3630.0 ±50.3 
638.2 ±38.4 
137.9 ±30.9 


126.7 ±25.6 

32.5 ± 4.0 
347.9 ± 128.6 
2200.0 ±98.2 
534.0 ±51.3 

86.6 ± 23.3 


1 192.4 ± 87.9 
766.5 + 52.1 
2563.6 + 711.7 
9140.0 + 324.7 
5235.6 ± 89.6 
946.1 +49.4 


888.6 ± 113.0 
1871.8 ± 186.5 

15030.4 ±335- 
1259.8 ± 289.7 


System 
VCC 2062 


FUV 

47.2 + 4.4 


NUV 
64.9 + 6.0 


163.0+ 114.7 


B 

234.5 + 39.8 


g 

251.3 +47.6 


V 

259.4 + 27.1 


268.6 + 63.0 


k 

281.4 + 38.7 


253.3 + 74.2 


271.7 + 215.0 


3.6 

210.0 + 52.0 


4.5 

185.0 + 55.0 


y.o 

1597.0 + 387.0 



TABLE 2 
Spectroscopic characteristics 



System 


Oxygen abundance 


EW(H£) 


A v 




12 + log O/H 


A 




Stephan's Quintet 2 


8.8 






Stephan's Quintet 5° 


8.7 


137-454 


1.3-2.3 


Arp 105S* 


8.4 


111 


0.8 


Arp 245N C 


8.6-8.7 


6-13 


2.0 


NGC 529 lN d 


8.4 


140 


<0.1 


NGC 7252NW C 


8.6 


30 


1.1 


VCC 2062^ 


8.6-8.7 


70 


0.1 



Note. — a=Duc et al. (in preparation), b=Duc & Mirabel 

(1994) , c=Duc et al. (2000), d=Duc & Mirabel (1998), e=Duc 

(1995) ,f=Ducet al. (2007). 




0.1 1.0 10.0 1 10 100 1000 

Wavelength [urn] Time [Myr] 

Fig. 3. — Modeling of region 2 in Stephan's Quintet. We use a mass ratio of young to old stars logr = -1. The fit yields no attenuation for the evolved 
population of galactic origin, and a timescale r = 37 X 10 6 years. The fit has been performed taking into account the FUV, NUV, B, V and R-band data. Given 
their uncertainties, the Ho- flux and the H/3 equivalent width have not been used to constrain the fits, but are included in the plot as an indication. Left: the spectral 
energy distribution of the best fit. The dotted line represents the spectral energy distribution of the population originating from the parent galaxies, the dashed 
line represents the population formed by the collision debris and the solid line represents the sum of both populations. The star symbols indicate the flux in the 
filters when convolved with the spectral energy distribution, taking into account the nebular emissions, whereas the red lines crossing the error bars represent the 
observations. Right: thex 2 as a function of the age and the attenuation of the young stellar component. The dark blue color indicates low^ 2 values and red high 
ones. The white contours indicate the value of the H/3 equivalent width in A, and the black ones the ratio of the Ha from the model to the observed value. The 
attenuation is taken into account. Finally the brown contours follow the \~ color scheme and give the probability that given model parameters can reproduce the 
observations. The cross shows which model SED has the highest probability to reproduce the observations. 
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Stephan's Quintet §2 Stephan's Quintet §2 




0.1 1 10 100 0.1 1 10 

Wavelength [jLim] Wavelength [/im] 



Fig. 5. — Left: Fitted SED of region 2 of Stephan's Quintet. Red: dust emission; blue: stellar spectrum without attenuation; cobalt: attenuated stellar spectrum; 
gray: total; green dots: model integrated within the broadbands. The different curves correspond to different values of Go: 1, 10 and 100. Right: Observed SED 
of region 2 of Stephan's Quintet (blue circles with error bars), compared to SED templates of nearby galaxies modeled by Galliano et al. (2008a) (red line). 
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Stephan's Quintet §5 
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Fig. 6. — Same as Figure 5 for Stephan's Quintet region 5. The value of Go is 

do not pretend that the far-IR SED will look as plotted. This 
fit simply shows the highest value of Go (here Go - 600) 
dominating the IR SED. It is equivalent to the parameter U max 
in the Dale et al. (2001) prescription. 

4.2. Arp 105S 

Arp 105S has the color and morphology of a blue compact 
dwarf galaxy. Due & Mirabel (1994) have shown that it is in 
fact a Tidal Dwarf Galaxy exhibiting prominent star-forming 
regions. The measured H/3 equivalent width of 1 1 1 A suggests 
a fairly young burst. 

Arp 105S is located towards the stellar halo of an elliptical 
galaxy (Fig. 1 of Boquien et al. 2009). So as to minimize the 
pollution by this galaxy, we have subtracted its profile using 
the ellipse and bmodel procedures in iraf. 

4.2.1. Stellar populations 

The best fit (presented in Figure 7, with a probability of 
67% to reproduce the observations) corresponds to a burst 
strength of log r — 2, an age of t — 4 x 10 6 years, a very short 
time scale of r = 10 6 years (i.e. an instantaneous burst) and 
an attenuation of 1.6 magnitudes for the population formed in 
the collision debris. The Ha flux is however overestimated 
and the attenuation is larger than the one determined by Due 
& Mirabel (1994). The ultraviolet flux is not perfectly ad- 
justed by the model, which may be due to a slight remaining 
contamination by the elliptical galaxy. 

To test whether the presence of an old population is compat- 
ible with the observations, we have also fitted the SED setting 
log r = — 1. The best fit (shown in Figure 8) then corresponds 
to a probability that the model can reproduce the observations 
of 63%. 

The parameters are close to the ones obtained for log r = 
2: t — 4 x 10 6 years, r = 10 6 years and a slightly lower 
attenuation. Values well below log r — - 1 are excluded by 
our models. 

4.2.2. Dust 

In Figure 9 we present the fit of Arp 105S with its best 
template and a dust model. 

The best template for this blue star-forming region, simi- 
larly to the region 5 of Stephan's Quintet, is the starburst BCD 
Mrk 930. The mismatch in the mid-IR tells about a relatively 
high dust content. 
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The rescaled IRS spectrum allows us to accurately constrain 
the PAH-to-dust mass ratio. However, the fit of the spectrum 
is not perfect. In particular, the 6.2 to 7.7 pm feature ratio 
is lower than the modeled one. This peculiarity could be ex- 
plained by a destruction of the smallest PAH in this partic- 
ularly active star-forming object (Fig. 14 of Galliano et al. 
2008b). 

4.3. Arp 245N 

As shown by Due et al. (2000), Arp 245N is a Tidal Dwarf 
Galaxy candidate located at the tip of a prominent stellar tidal 
tail (Fig. 2 of Boquien et al. 2009). The optical images of 
Arp 245N show without any possible doubt the presence of 
an evolved population originating from the parent galaxies. 
The TDG exhibits several bluer regions, with Ha, UV and 
mid-IR counterparts, indicating the presence of on-going star 
formation. 

4.3.1. Stellar populations 

The best fit (shown in Figure 10) corresponds to a burst 
strength of only logr = -0.75, t — 115 x 10 6 years, a short 
timescale of t = 10 6 years, an attenuation of 1.1 magnitude 
for the evolved population and 1 .5 magnitudes for the popu- 
lation formed in the collision debris. 

Only a narrow range of parameters has a non negligible 
probability to reproduce the observations. 

Given the strong constraints given by the fit, we could for 
this system analyze the effect of changing the initial mass 
function. The best fit using a Scalo (1986) IMF is shown in 
figure 11. 

The burst strength remains unchanged log r = -0.75 and 
the other parameters change slightly: t = 100 x 10 6 years; t = 
10 6 years; attenuations of 0.6 magnitude and 1.2 magnitudes 
for the old and young stellar populations. 

The small offset towards younger ages is due to the fact that 
the Scalo (1986) IMF creates less massive stars, which have 
the strongest effect in the ultraviolet. 

We note that with the Scalo (1986) initial mass function, 
the IRAC bands at 3.6 pm and 4.5 pm are better reproduced 
than with the Salpeter (1955) initial mass function. With the 
Scalo (1986) IMF more intermediate and less very high or 
very low mass stars are formed than with the Salpeter (1955) 
IMF which leads to lower 3.6/optical and 4.5/optical ratio. 

The spectral energy distribution is best fitted assuming a 
Scalo (1986) initial mass function. If we assume that it is 
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Fig. 8. — Modeling of Arp 105S with a mass ratio log r = — 1, an attenuation of 1.4 magnitudes for the evolved component and 1.5 magnitudes for the population 
formed in collision debris and a timescale r = 10 6 years. The fit has been performed on the FUV, NUV, B, V, R, J, H and K bands. See Figure 3 for additional 
details. 
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Fig. 11. — Modeling of Arp 245N with a Scalo (1986) initial mass function, a mass ratio logr = -0.75, an attenuation of 0.6 magnitudes for the evolved 
component and a timescale r = 1 X 10 6 years. The fit has been performed on the FUV, NUV, B, V, R, J, H and K bands. See Figure 3 for additional details. 



the correct initial mass function, it means that the 3.6 fim and 
4.5 fim bands are largely dominated by the stellar flux and 
contain very little pollution by a dust component. Unsurpris- 
ingly, Boquien et al. (2009) also found that the very small 
grains emission was small in Arp 245N compared to the PAH 
emission. 

4.3.2. Dust 

In Figure 12 we present the fit of Arp 245N with its best 
template and a dust model. 

The best template in the UV-to-IR regime corresponds to 
NGC 520, a starburst merger with composite stellar popula- 
tions, which is consistent with our modeling of its stars. In- 
terestingly, Arp 245N shows weaker 6.2 to 7.7 fim and PAH 
emission compared to the one of NGC 520. 

4.4. NGC 5291N 

NGC 529 IN is the most luminous star forming region in 
the collisional ring around NGC 5291 (Fig. 3 of Boquien 
et al. 2009). Previous works by Due & Mirabel (1998); Hig- 
don et al. (2006); Boquien et al. (2007) suggest that the star 
forming regions in this system are probably very young and 
dynamical arguments hint that there might not be any under- 
lying evolved population. The nebular lines are particularly 
strong in this system as has been shown by Due & Mirabel 
(1998), which makes a very accurate SED modeling more dif- 
ficult. 

4.4.1. Stellar populations 

The best fit (shown in Figure 1 3) indicates a rather young 
burst, t — 5 x 10 6 years with a short timescale t — 10 6 years, 
an attenuation of 1.1 magnitudes and a high burst strength 
log r — 2. The flux in the J band is slightly overevaluated 
by the model. This may be due to the uncertainties on the 



emission line fluxes. For a young starburst this band is dom- 
inated by the Pa-/? line at 1.28181 fim and by the [Ne II] line 
at 1.28 fim. 

At wavelengths above 1.3 fim, the observed fluxes are 1.2- 
2.5 times higher than those predicted by the model. In prin- 
ciple, this could trace the presence of an old stellar popu- 
lation emitting in the near-infrared. To check this, we set 
log r = -1. As shown in Figure 14, if the fit is better for the K 
band, it is not improved in the other near-infrared bands and it 
is significantly degraded in the optical bands, the fluxes being 
systematically under-evaluated. If not due to stellar contin- 
uum, the excess of light in the near-infrared regime should 
already be attributed to dust emission, usually considered to 
become dominant at longer wavelength. This hypothesis will 
further be discussed in section 5.2.3. 

4.4.2. Dust 

In Figure 15 we present the fit of the SED of NGC 529 IN 
with its best template and a dust model. 

The best galaxy template in the UV-to-optical regime is the 
one of UM 448 (Mrk 1304), a peculiar starburst galaxy with 
a relatively low metallicity (12+log(0/H)=8) but which ex- 
hibits faint PAH emission (Engelbracht et al. 2005). Like for 
the other active star forming regions in our sample, the emis- 
sion in the mid-IR is much stronger. NGC 529 IN has more 
prominent PAH features. 

We fit its SED with a dust model corresponding to three 
different cases: 1. Ag = 0.2 as observed by Due & Mirabel 
(1998) which corresponds to the lower infrared emission, 2. 
with a SMC attenuation curve which corresponds to the mid- 
dle infrared emission and finally 3. using the full absorbed 
power constraint. The latter two cases fit similarly the data, 
however Ab = 0.2 is too low to reproduce the 24 fim data. 

4.5. NGC 7252NW 
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Fig. 12. — Same as Figure 5 for Arp 245N. The black line represents the Spitzer rescaled IRS spectrum (Boquien et al. 2009). The values of Go are 1, 10 and 
100. 
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Fig. 13. — Modeling of NGC 529 IN with a mass ratio log r = 2, an attenuation of 1.1 magnitudes for the population formed in collision debris and a timescale 
r = 1 X 10 6 years. The fit has been performed on the FUV, NUV, B, V, R, J, H and K bands. See Figure 3 for additional details. 



NGC 7252NW is a star-forming region located at the tip of 
a stellar tidal tail emanating from an old merger remnant (Fig. 
4 in Boquien et al. 2009). 

4.5.1. Stellar populations 

The best fit shown in Figure 16 reproduces the observa- 
tions with a probability of 73% with the following parame- 
ters: a rather small burst strength of log r = -0.75, an age 
t = 300 x 10 6 years, a timescale r = 175 x 10 6 years an at- 
tenuation of 0.9 magnitude for the evolved population and an 
attenuation of 0.2 magnitude for the young population, which 
is significantly smaller than what was observed spectroscop- 
ically. In all cases, reasonable fits yield an attenuation lower 
than 1 magnitude. 

There is a wide range of possible timescales that can re- 
produce the observations with a probability higher than 50% 



as we can see in Figure 17. A nearly instantaneous starburst 
leads to an age a bit younger than 60 x 10 6 years. On the con- 
trary, for a nearly constant star formation, the best fit occurs 
at an age of about t = 750 x 10 6 years. However this is an 
upper limit at fixed r and attenuation. 

As the merger is old, a long timescale is not unexpected. 
Though, the attenuation is questionably small whereas a 
strong excess can be seen in the infrared. 

4.5.2. Dust 

In Figure 18 we present the fit the SED of NGC 7252NW 
with its best template and a dust model. 

The best template corresponds to NGC 5253, a nearby 
dwarf galaxy which has exhibited starbursting episodes dur- 
ing the last few hundred million years, and where faint PAH 
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Fig. 14. — Modeling of NGC 529 IN with a mass ratio log r = -1, an attenuation of 0.9 magnitude for the population formed in collision debris, no attenuation 
for the evolved population and a timescale r = 1 X 10 6 years. The fit has been performed on the FUV, NUV, B, V, R, J, H and K bands. See Figure 3 for additional 
details. 
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Fig. 15. — Left: Fitted SED of NGC 5291. The different models correspond to (1) the case where we fix Ag = 0.2 as observed spectroscopically by Due & 
Mirabel (1998) (lowest IR emission) with Go = 10 5 and a starburst attenuation law, (2) we adopt the SMC attenuation curve (middle IR emission) with Go = 445, 
(3) where we use the full absorbed power constraint (highest IR emission) with Go = 77. The attenuation in cases (2) and (3) is derived from the mode. Right: 
Same as Figure 5 for NGC 5291. The black line represents the rescaled Spitzer IRS spectrum (Higdon et al. 2006). 



features were detected (Beirao et al. 2006). Whereas the emis- 
sion at 8 yum of NGC 7252NW is much higher, probably due 
to the presence of stronger PAH, its 24 /jm is surprisingly 
lower. The dust model also shows inconsistent results. The 
attenuation measured from the line ratio is so low that most of 
the IR emission would be coming from PAH. This is highly 
questionable and significantly smaller than what was observed 
spectroscopically. This will be discussed in section 5.2.2. 

4.6. VCC2062 

VCC 2062 has been described as an old TDG candidate, as 
its parent galaxy has now evolved into a relaxed early-type 
galaxy (Due et al. 2007). It is linked to the latter by a star-less 



HI bridge. 

Unfortunately this system does not have any near-infrared 
image available. As no diffuse stellar component that would 
hint about the presence of an evolved stellar population can be 
seen, we assume log r = 2. Note that a preliminary spectral 
energy distribution fit was published in Due et al. (2007). 

4.6.1. Stellar populations 

In Figure 19 we present the best fit of the spectral energy 
distribution. In agreement with the spectroscopic observa- 
tions, the attenuation is weak, typically a few tenths of a mag- 
nitude. The best fit yields an age of t — 1 15 x 10 6 years and a 
nearly instantaneous burst: t = 10 6 years. 
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Fig. 16. — Modeling of NGC 7252NW with a mass ratio log r = -0.75, 0.9 magnitude of attenuation for the evolved population and a timescale r = 175 X 10 6 
years. The fit has been performed on the FUV, NUV, B, V, R and N2 bands. See Figure 3 for additional details. 
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Fig. 17. — Modeling of NGC 7252NW with a mass ratio log r = -0.75, 0.9 magnitude of attenuation for the evolved population and an attenuation of 0.2 
magnitude for the population formed in the collision debris. The fit has been performed on the FUV, NUV, B, V, R and N2 bands. See Figure 3 for additional 
details. 
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To verify whether this result is still valid in the subregions 
of VCC 2062 where no more star formation can be detected, 
we have also performed a similar fit on the North-East part of 
the galaxy; there no Ha emission can be detected. Qualita- 
tively the results are very similar, with an age t = 270 x 10 6 
years and an attenuation of 0.3 magnitude. The timescale is 
also very short: t = 2x 10 6 years. 

We see in Figure 20 that whatever the timescale is, the age 
of the burst cannot be younger than approximately t = 100 x 
10 6 years. This age changes little for a timescale shorter than 
a few tens of million years. However in this configuration the 
Ho- emission is too weak. To obtain a flux similar to the one 
observed, a longer timescale is necessary, about t = 50 x 10 6 
years and an age of about t = 200 x 10 6 years. Beyond, the 
Ho- flux becomes significantly too high even though the fits 
remain proper. 

Because of the uncertainties on the fluxes and the lack of 
near-infrared data beyond the z' band - which flux is particu- 
larly uncertain -, it is hard to have definitive conclusions re- 
garding this star forming region. The fits we have performed 
indicate that star formation started between t = 200 x 10 6 and 
t = 300 x 10 6 years ago. The timescale is short compared to 
the age, at most a few tenths of it. 

4.6.2. Dust 

In Figure 21 we present the fit of the SED of VCC 2062 
with its best template and a dust model. 

Like for SQ-2, the best galaxy template of VCC 2062 corre- 
sponds to the spiral M5 1 . Similarly, the fit is reasonably good 
from the UV to the mid-IR. 

5. DISCUSSION 
5.1. Stellar populations in collision debris 
5.1.1. Old stellar populations 

The stellar populations in tidal debris consist of an old stel- 
lar population expelled from the parent galaxies and a young 
one born in situ. The mass fraction of each type of stars had 
so far been determined qualitatively, based on the gas-richness 
and colors of the tidal tails: the more gas rich and the bluer 
they are, the higher the apparent contribution of young stars. 
First attempts to determine the stellar populations in a more 
quantitative way were made by Fritze-v. Alvensleben & Due 
(1998) based on B - V, V - R and V - K colors. In this paper, 
we have used a much wider range, from the UV to the mid-IR 
to get more precise insights. Our study reveals a variety of 
properties in those objects. 

Among the selected systems, four are compatible with 
a population constituted exclusively in the collision debris, 
without any underlying evolved stellar population. The con- 
clusion is drawn from the best fit but it is also based on a series 
of tests as, for instance, increasing log r as was shown in sec- 
tion 4.2. In NGC 5291N, our result is corroborated by numer- 
ical simulations (Bournaud et al. 2007) indicating that there 
should not be any population of galactic origin located in the 
collisional ring. Regarding Arp 105S, we cannot completely 
exclude the possibility of the presence of an evolved popula- 
tion even though it is unlikely. The presence of the nearby 
elliptical galaxy makes accurate photometric measurements 
difficult. In the case of SQ-5 and VCC 2062, even though fits 
lean towards the absence of an older stellar population, the 
availability of near-infrared observations is crucial to draw a 
more definitive answer. 



Other systems such as SQ-2, Arp 245N and NGC 7252NW 
most likely contain an evolved stellar population, as was ex- 
pected. In the case of Arp 245N, it was clearly detected and 
constrained thanks to near-infrared observations in the J, H 
and K bands. 

For young objects, a ratio as low as 10% is sufficient for 
the young stellar population to dominate the spectral energy 
distribution, masking the presence of an underlying evolved 
stellar population as can be seen in Figure 22. 

This 10% limit depends of course strongly on the star for- 
mation history and on the quality of the available observa- 
tions to disentangle the emission from both components. For 
a mass ratio log r = -1 at 5 x 10 6 years, the young population 
accounts for at least 75% of the flux, the minimum being in 
the H band. Yet, for the systems in which there is clearly an 
evolved stellar population, its emission in the near-infrared is 
similar to the error bars for the currently available data. Thus, 
for a given age, the detectability limit depends on the error 
bars in the near-infrared. In the present case, for a 5 x 10 6 
years old burst, it is about 10%. 

5.1.2. Recent star formation 

The recent, in-situ, star formation history in collisional de- 
bris may be be inferred comparing various tracers of star for- 
mation: ultraviolet, Ha or mid-/far-infrared. 

Each tracer has its own characteristics. Ultraviolet is sen- 
sitive to the photospheric emission of massive stars over a 
timescale of a few 100 x 10 6 years, Ha to the ionized gas 
surrounded star forming regions over a timescale of about 
10 x 10 6 years and finally infrared is a result of the radiation 
of dust absorbing emission from massive young stars. 

The availability of several bands permits not only to con- 
strain the recent star formation history but also to obtain a 
more precise estimate of the star formation rates. However 
they are still sensitive to some degeneracies, depending on 
the observations available, and thus the fits must be cautiously 
interpreted. We compare the SFR obtained from the SED fit- 
ting, presented in Table 3 to the ones published in Boquien 
et al. (2009) which use the conversion factors published by 
Kennicutt (1998) for the UV and the Ha emission and by Wu 
et al. (2005) for the 8 pm. The SFR given by the models are 
the "true" instantaneous SFR. The ones published in Boquien 
et al. (2009) assumed a given star formation history. 

We notice that the model strongly underestimates the SFR 
from the observations for VCC 2062 if we assume a sin- 
gle recent burst. This is most likely due to several factors: 
the lack of observations in near-infrared beyond the z' band 
and a star formation history that is not properly reproduced 
by our set of models. An older age of about 300 Myr with 
a timescale around 100 Myr, which is not ruled out by the 
model is probably more realistic as it would take into ac- 
count more accurately a long term low level star formation. 
For the other models, except for NGC 529 IN, SFR(FUV) 
systematically underestimates SFR(model). In the case of 
SQ-5, SFR(model)/SFR(FUV) even reaches ~ 32. This is 
chiefly an attenuation effect: spectroscopic observations of 
NGC 529 IN show that the attenuation is very small (Due & 
Mirabel 1998) and SQ-5 is particularly bright at 8 pm com- 
pared to FUV and has a relatively high attenuation. Regard- 
ing ionized gas, SFR(model)/SFR(Ha) is comprised between 
0.73 for NGC 7252 and 2.8 for SQ-5. No trend can be seen 
with the star formation history or the attenuation hinting at a 
complex combination of several factors. Finally, SFR(model) 
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Fig. 20. — Modeling of VCC 2062 with a mass ratio log r = 2 and an attenuation of 0.4 magnitude for the population formed in the collision debris. The fit has 
been performed on the FUV, NUV, u', B, g', V, r', R, i' and z' bands. See Figure 3 for additional details. 
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Fig. 21. — Same as Figure 5 for VCC 2062. The values of Go are 1,10 and 100. 




and SFR(8.0) - the star formation rate obtained from the PAH 
emission in the Spitzer 8.0 fim band - provide similar re- 
sults with two exceptions: SQ-5 reaches a ratio of ~ 4.5 and 
Arp 245N presents a strong mid-infrared emission with a ratio 
of 0.28. 

Even though the standard SFR estimators often yield results 
that have the correct order of magnitude, an error of several 
times is common with extreme cases of up to several tens. 
This shows that as SED fits allow to model the star formation 
history and the attenuation they permit to obtain more accu- 
rate results. The SFR from the Ha provides relatively good 
estimates of the real star formation rate, with a maximum ex- 
cursion of 2.44. This is not surprising as the timescale of Hot 
emission is small and therefore more closely traces recent star 
formation. 



5.2. Dust in collision debris 

As star-formation starts in dust cocoons, the dust compo- 
nent should be an important ingredient of collisional debris 
and indeed it has been observed in this environment either in 
absorption, as dust lanes, that can be seen for instance in the 
region B of Stephan's Quintet (Lisenfeld et al. 2004), or di- 
rectly in emission (Higdon et al. 2006; Boquien et al. 2009). 
Various types of dust may be probed by our data, as discussed 
below. 

5.2.1. A qualitative comparison with UV-to-IR galaxy SED 
templates 

In section 4, we have shown the comparison of the observed 
SED of each star forming region with the SED of nearby 
galaxies modeled by Galliano et al. (2008a). Their 35 tem- 
plates were scaled to match the UV-to-near-IR (up to 4.5 //m) 
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Fig. 22. — Evolution of the contribution of the flux in the young population 
relatively to the total flux in each band used during this study. The burst is 
nearly instantaneous et there is no attenuation. The solid line represents an 
age of 5 X 10 6 years, the mixed one an age of 50 X 10 6 years and the dashed 
one an age of 550 X 10 6 years. The bands are FUV, V, H, 3.6 fim, and 8.0 fim. 



observed SED of the galaxies presented in this study. We kept 
the template minimizing the^ 2 . 

For Arp 105S, Stephan's Quintet 5, NGC 7252NW and 
NGC 529 IN, the templates that match best the stellar spec- 
trum are templates of low-metallicity, galaxies. Indeed, the 
stellar SED of these objects is very blue. However, when 
comparing the fitted template to the mid-IR observations, the 
observations present an excess emission. This indicates that 
the PAH abundance is higher than in the templates. High PAH 
abundance is the sign of a metal rich ISM. Therefore, those 4 
intergalactic star forming regions have a young stellar popu- 
lation typical of a blue compact galaxy, but have the mid-IR 
emission of a solar metallicity galaxy. 

For the rest of the sample, Arp 245S, Stephan's Quintet 2, 
and VCC 2062, are relatively well matched by SED templates 
of spiral galaxies, because their stellar spectrum is redder. 

5.2.2. The PAH-to-dust mass ratio and the metallicity 

The PAH-to-dust mass ratio is a function of ISM metallic- 
ity. Madden et al. (2006) and Wu et al. (2006) have shown 
that the PAH feature strength to continuum intensity ratio 
is correlated both with the metallicity of the ISM and the 
[Ne m]/[Ne n] line ratio. The latter is a tracer of very young 
stellar populations, and therefore an indicator of the hardness 
of radiation field. Madden et al. (2006) proposed that the 
paucity of PAH in low-metallicity environments was a conse- 
quence of their large-scale destruction by the hard penetrating 
radiation. Indeed, those environments contain more ionizing 
radiation, and their ISM has a lower opacity, due to the lower 
dust abundance. 



Galliano et al. (2008a) showed that although PAH destruc- 
tion is important, the trend is governed by the delayed injec- 
tion of PAH by their progenitors, the AGB (asymptotic giant 
branch) stars. AGB stars are believed to be the main source of 
PAH. AGB stars begin to inject their elements and dust after 
their lifetime of ^ 400 Myr, while massive stars inject their 
elements and dust (non-PAH carbon and silicate dust) into the 
ISM after only a few 10 6 years. The metallicity being a func- 
tion of the age of the galaxy, the delayed injection of PAH 
by AGB stars explains the deficit of PAH in low-metallicity 
environments. 

Although the origin of this trend is still controversial, the 
trend in itself is consensual. We can therefore assume that the 
PAH-to-dust mass ratio scales with the metallicity of each of 
our intergalactic star forming regions. 

As a test of the attenuation obtained with the photo- 
spheric +nebular model, Figure 23 compares the values of the 
PAH-to-total dust mass (comprising the PAH, the Very Small 
Grains and the Big Grains) ratio /pah to the range of values 
consistent with the observed trend with metallicity (Galliano 
et al. 2008a). We emphasize that this value of /pah relies on 
the accuracy of the attenuation estimate. Here, we do not con- 
sider this determination of f PA H as an actual measure of the 
PAH-to-dust mass ratio, but we rather compare this derived 
fpAH to its expected value, taking the trend of PAH abundance 
with metallicity for granted. From this comparison, there are 
three configurations. First, if these two values are in agree- 
ment, it provides an additional indication that the attenuation 
estimate, as well as the fitted stellar populations, might be ac- 
curate. Second, if the derived /pah is much higher than its ex- 
pected value, it is the indication that the attenuation estimate 
implies a gross underestimate of the IR luminosity. It there- 
fore suggests that regions contributing to the IR luminosity 
were not probed in attenuation by the photospheric+nebular 
fit. These regions are likely embedded star forming regions, 
contributing to the IR power, but optically thick to visible re- 
combination lines. Finally, if the derived value of /pah was 
lower than its expected value, it would mean that we overesti- 
mated the attenuation. This situation would be more difficult 
to justify from a physical point of view. However, we never 
encounter this situation in our sample. 

Most objects are within the expected range, except 
NGC 529 IN and NGC 7252NW, in the case where we adopt 
Ab - 0.2 as was observed spectroscopically by Due & 
Mirabel (1998). The fact that these two objects are outliers 
is an indication that the attenuation estimate does not account 
for enough IR emission. It is obvious from their SED them- 
selves (Figures 15 and 18, left panel), where we see that for 
those two cases, the IR emission comes essentially from the 
PAH. This discrepancy allows us to estimate that there is a 
missing IR luminosity, likely coming from embedded star 
forming regions. For NGC 5291N and NGC 7252NW, this 
missing IR luminosity is about an order of magnitude higher 
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Fig. 23. — Comparison the derived PAH-to-dust mass ratio to their theoreti- 
cal value. The points with error bars are the PAH-to-dust mass ratios derived 
from the SED fit, as a function metallicity. The gray curve is the theoretical 
value of this ratio, from chemical evolution considerations (Galliano et al. 
2008a). 

than the power derived from the attenuation estimate. Table 4 
summarizes those parameters. 

We also note that there is an anticorrelation between 
M(PAH)/M(dust) and Av. A possibility would be that it is 
due to the color of the ISRF. However, it is lower than the 
variations we are probing. For example, Fig. 3 of Galliano 
et al. (2008a) demonstrates the effect of the color of the ISRF 
on the PAH spectrum, for two extreme cases. The varia- 
tion is at most a factor of 3, and it is certainly lower in this 
sample: we are probing variations of 2 order of magnitude 
in M(PAH)/M(dust). The reason of this anticorrelation most 
likely reflects the systematic effect of the determination of 
Av. If we consider the V band flux over 8 fim in our sam- 
ple, it does not vary drastically: F Y (V)/F V (8) 0.1, SQ5 
being lower (0.02) and NGC 245 higher (0.25). At first or- 
der, the value of Av will determine the fraction of absorbed 
luminosity in the UV-visible, relative to the V band. Since, 
the V band to 8 fim flux does not vary a lot, Av will deter- 
mine the absorbed luminosity relative to the 8 fim band. Fi- 
nally, the value of M(PAH)/M(dust) will depend directly on 
the 8 fim to absorbed luminosity ratio. Therefore, when Av is 
low, 8 yizm/absorbed luminosity is going to be high, and there- 
fore, M(PAH)/M(dust) is going to be high. Fig. 23 is a way to 
help discriminate the values of Av which are clearly too low 
compared to the known 8 fim emission. A very low value of 
Av (e.g. NGC 5291 Av=0.2) would require all the observed 
luminosity to be put in the PAHs. There would be no room 
for far-IR dust emission (Fig. 14), which is difficult to under- 
stand theoretically, and has never been observed at the scale 
of a galaxy. 

5.2.3. Dust contribution to the 3.6 fim and 4.5 fim bands 

The modeling of the fluxes at 3.6 fim and 4.5 /im by the 
photospheric and nebular models is not always very satisfy- 
ing, in particular for the regions for which the luminosity is 
dominated by the population formed in the collision debris. 
To exhibit the physical cause, we have plotted the ratio of 
the modeled to the observed bands versus the observed 8 fim 
band emission (which traces star formation) normalized to the 
emission of the model at 3.6 fim, which can be seen as a spe- 
cific star formation rate. In Figure 24 we see that as the spe- 
cific star formation rate increases, the fit of the 3.6 fim and 



4.5 fim bands degrades progressively. The Pearson correla- 
tion coefficient is r — -0.71 (resp. r = -0.72) for the 3.6 fim 
(resp. 4.5 fim) band. 

Smith & Hancock (2009) have made an extensive study on 
the origin of the non-stellar 4.5 fim excess in dwarf galaxies 
investigating several possibilities. They have shown this ex- 
cess is likely a combination of nebular emission (Br-c line 
and continuum) and high optical depth - which would red- 
den the 3.6 fim to 4.5 fim color -, making the presence of 
dust emission unnecessary. However we cannot rule out that 
the contamination of the 3.6 fim and 4.5 fim bands by dust 
emission is a possibility in star forming regions in collision 
debris. A high optical depth may be a factor in NGC 5291N 
and NGC 7252NW as we have shown that there are hints 
of embedded star forming regions. Our UV to near-infrared 
model already takes into account the nebular emission, both 
the Br-ff line and the continuum, so the detected excess is 
unlikely to be of nebular origin in its entirety. This leads us 
to interpret this excess as the effect of the emission of the 
dust (the continuum as well as the emission lines) excited 
and heated by massive stars. This emission is probably com- 
ing from very small grains fluctuating in temperature, rather 
than hot grains, because the equilibrium temperature of the 
latter (^700K) would correspond to grains too close to stel- 
lar sources (^1 a.u. of an 04 star). When the specific star 
formation rate is high, there are proportionally many more 
young stars that are able to make the small grains fluctuate 
to higher temperatures than the rest of the stellar population. 
As a consequence, the dust emission at short wavelength be- 
comes non-negligible compared to the stellar emission. As a 
consequence, the 3.6 fim and 4.5 fim bands cannot be good 
tracers of the mass of the evolved stellar population, espe- 
cially in regions that have a high specific star formation rate. 
In addition, the objects that have the highest excess at 3.6 fim 
and 4.5 fim are also those with a negligible evolved stellar 
population that can contribute to those bands, making it eas- 
ier for the dust to dominate. An uncertainty may be intro- 
duced by the r parameter. It should have little influence as the 
J, H and K near-infrared bands strongly constrain the mod- 
eled SED in longward bands in wich the stellar emission is 
dominated by the same populations. However, when no near- 
infrared band is available shortward of 3.6 fim, the r parameter 
is much harder to constrain and the uncertainty in the 3.6 fim 
and 4.5 fim bands is much higher. 

6. SUMMARY AND CONCLUSIONS 

In this study we have compiled and modeled the spectral 
energy distribution of a sample of 7 star forming regions in 
collision debris. We have carried out this study in three steps. 
First, we have adjusted the SED from the ultraviolet to the 
near-infrared (if available, otherwise ultraviolet to optical) 
with a photospheric +nebular model, yielding information on 
the star formation rate, attenuation, and star formation his- 
tory. In a second step, we have fitted their UV to near-IR SED 
with templates of well known galaxies and then compared the 
mid-infrared part of their SED to infer information on the dust 
content. Finally, a dust model was made to obtain estimates of 
the dust-to-PAH mass ratios. We have obtained the following 
main results: 

1 . We confirm that several star forming regions seem com- 
patible with a complete absence of any evolved stel- 
lar population: Stephan's Quintet 5, Arp 105S and 
NGC 529 IN. This makes them ideal to study star for- 
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TABLE 4 

Parameters relative to the mid-IR SED fit. 



Object 


M PAH /M dust 


Absorbed power [L ] 


Missing L IR [L ] 


Arp 105S 


0.003 ± 0.001 


3.1 x 10 10 




Arp 245N 


0.028 ± 0.001 


10 9 




Stephan's Quintet #2 


0.015 ± 0.001 


10 9 




Stephan's Quintet #5 


0.030 ± 0.002 


5.1 x 10 8 




NGC 7252NW 


0.21 ±0.01 


10 8 


10 8 -2x 10 9 


VCC2062 


0.10 ±0.02 


10 7 




NGC5291N(A B = 0.2) 


0.78 ± 0.01 


2.9 x 10 8 


^ 2 x 10 9 - 8 x 10 10 


NGC 529 IN 


0.0046 ± 0.0001 


10 10 




NGC 529 IN (SMC ext.) 


0.023 ± 0.0004 


3.0 x 10 9 




Galaxy 
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Fig. 24. — Left: Ratio of 3.6 /im model flux to the 3.6 /im observed flux versus the specific star formation rate. Right: same as left but for the 4.5 /im band. The 
fluxes are obtained from photospheric and nebular only models, not from the dust models. 



mation isolated from the influence of any evolved stel- 
lar population. 

2. The UV-to-optical SED of these star forming regions 
exhibiting high burst strengths are best fit with tem- 
plates of starbursting dwarf galaxies, while their near 
to mid-IR SED resemble more that of dusty, metal- 
rich, star-forming regions in spiral disks. The colli- 
sion debris dominated by an old stellar component, like 
Stephan's Quintet 2, show a good fit with spiral disks 
from the UV to the mid-infrared. 

3. The 3.6 yum and 4.5 fan bands have a luminosity that 
cannot be explained by the stellar luminosity alone. 
This excess is correlated with the star formation in- 
tensity. The contamination most likely originates from 
very small grains. As such, those bands should not be 
used to trace the stellar mass in collision debris unless 
the star formation intensity is weak. 



4. The 8 //m emission provides an additional consistency 
check, showing that most of our attenuation estimates 
and subsequent stellar population synthesis are reliable. 
We find discrepancies in two cases: NGC 5291N and 
NGC 7252NW. The nature of these discrepancies indi- 
cates that there is a significant fraction of the infrared 
power, probably originating in embedded star form- 
ing regions, that was not probed through the photo- 
spheric+nebular fit, in those sources. 

Deeper and additional data, especially in the near-infrared, 
are needed to constrain more accurately the mass of the var- 
ious stellar populations as the youngest population tends the 
dominate the radiation field even when representing a small 
fraction of the total stellar mass. A tentative constraint of the 
IMF has been obtained on one object. More accurate data are 
required to get significant results for the other regions in our 
sample. 

Furthermore, the interstellar radiation field intensity cannot 
be determined accurately without data in the far-infrared. This 
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makes a detailed study of the dust properties (such as the grain 
size distribution for instance) purely speculative. The recently 
launched Herschel Space Observatory is the key instrument to 
perform such a study, providing crucial far-infrared data while 
retaining a good spatial resolution to resolve star forming re- 
gions in collision debris from their nearby environment. 

We would like to thank the anonymous referee for the 
useful comments that have helped improve and clarify this 

REFE] 

Amram, P., Mendes de Oliveira, C, Plana, H., Balkowski, C, & Hernandez, 

O. 2007, A&A, 471,753 
Anders, P. & Fritze-v. Alvensleben, U. 2003, A&A, 401, 1063 
Appleton, P. N., Xu, K. C, Reach, W., et al. 2006, ApJ, 639, L51 
Beirao, P., Brandl, B. R., Devost, D., et al. 2006, ApJ, 643, LI 
Boquien, M., Due, P.-A., Braine, J., et al. 2007, A&A, 467, 93 
Boquien, M., Due, P.-A., Wu, Y, et al. 2009, AJ, 137, 4561 
Bournaud, R, Due, P.-A., Amram, P., Combes, R, & Gach, J.-L. 2004, 

A&A, 425, 813 

Bournaud, R, Due, P.-A., Brinks, E., et al. 2007, Science, 316, 1166 

Bournaud, R, Due, P.-A., & Masset, F. 2003, A&A, 411, L469 

Braine, J., Due, P.-A., Lisenfeld, U., et al. 2001, A&A, 378, 51 

Braine, J., Lisenfeld, U., Due, P.-A., & Leon, S. 2000, Nature, 403, 867 

Calzetti, D. 2001, PASP, 113, 1449 

Cluver, M., Appleton, P., Guillard, P., et al. 2010, ApJ 

Cortese, L., Gavazzi, G., Boselli, A., et al. 2006, A&A, 453, 847 

Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & Kolhatkar, S. 

2001, ApJ, 549,215 
de Grijs, R., Lee, J. T., Mora Herrera, M. C, Fritze-v. Alvensleben, U., & 

Anders, P. 2003, New Astronomy, 8, 155 
de Mello, D. R, Smith, L. J., Sabbi, E., et al. 2008, AJ, 135, 548 
Due, P.-A. 1995, PhD thesis, Universite Paris VI 
Due, P.-A., Bournaud, R, & Masset, F. 2004, A&A, 427, 803 
Due, P.-A., Braine, J., Lisenfeld, U., Brinks, E., & Boquien, M. 2007, A&A, 

475, 187 

Due, P.-A., Brinks, E., Springel, V., et al. 2000, AJ, 120, 1238 
Due, P.-A. & Mirabel, I. F. 1994, A&A, 289, 83 
Due, P.-A. & Mirabel, I. F. 1998, A&A, 333, 813 

Engelbracht, C. W., Gordon, K. D., Rieke, G. H., et al. 2005, ApJ, 628, L29 
Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Fritze-v. Alvensleben, U. & Due, P.-A. 1998, in The Magellanic Clouds and 
Other Dwarf Galaxies, Proceedings of the 

Bonn/Bochum-Graduiertenkolleg Workshop, held at the Physikzentrum 

Bad Honnef, Germany, January 19-22, 1998, Eds.: T. Richtler and J.M. 

Braun, Shaker Verlag, Aachen, ISBN 3-8265-4457-9, p. 141-145,, ed. 

T. Richtler & J. M. Braun, 141-145 
Galliano, R, Dwek, E., & Chanial, P. 2008a, ApJ, 672, 214 
Galliano, R, Madden, S. C, Tielens, A. G. G. M., Peeters, E., & Jones, A. P. 

2008b, ApJ, 679, 310 
Gerhard, O., Arnaboldi, M., Freeman, K. C, & Okamura, S. 2002, ApJ, 580, 

L121 

Gordon, K. D., Clayton, G. C, Misselt, K. A., Landolt, A. U., & Wolff, M. J. 

2003, ApJ, 594, 279 
Guhathakurta, P. & Draine, B. T. 1989, ApJ, 345, 230 
Hancock, M., Smith, B. J., Struck, C, et al. 2007, AJ, 133, 676 
Hancock, M., Smith, B. J., Struck, C, Giroux, M. L., & Hurlock, S. 2009, 

AJ, 137, 4643 



manuscript. 

UL acknowledges financial support by the Spanish Science 
Ministry under grant AYA 2007-67625-C02-02 and by the 
Junta de Andalucfa. VC would like to acknowledge partial 
support from the EU ToK grant 39965 and FP7-REGPOT 
206469. 

Facilities: Akari, Spitzer, GALEX, CAO:2.2m, CFHT, 
ESO:3.6m, KPNO:2.1m, NTT. 



Higdon, S. J., Higdon, J. L., & Marshall, J. 2006, ApJ, 640, 768 
Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 
Lisenfeld, U., Braine, J., Due, P.-A., et al. 2004, A&A, 426, 471 
Lisenfeld, U., Braine, J., Due, P.-A., et al. 2002, A&A, 394, 823 
Lopez-Sanchez, A. R., Esteban, C, & Rodriguez, M. 2004, ApJS, 153, 243 
Madden, S. C, Galliano, R, Jones, A. P., & Sauvage, M. 2006, A&A, 446, 
877 

Maraston, C. 2005, MNRAS, 362, 799 

Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 
Mendes de Oliveira, C, Cypriano, E. S., Sodre, Jr., L., & Balkowski, C. 

2004, ApJ, 605, L17 
Mendes de Oliveira, C. L., Temporin, S., Cypriano, E. S., et al. 2006, AJ, 

132,570 

Mezger, P. O. 1978, A&A, 70, 565 
Milgrom, M. 2007, ApJ, 667, L45 

Moles, M., Sulentic, J. W., & Marquez, I. 1997, ApJ, 485, L69+ 
Ryan-Weber, E. V., Meurer, G. R., Freeman, K. C, et al. 2004, AJ, 127, 
1431 

Sakai, S., Kennicutt, Jr., R. C, van der Hulst, J. M., & Moss, C. 2002, ApJ, 
578, 842 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11,1 
Schombert, J. M., Wallin, J. R, & Struck-Marcell, C. 1990, AJ, 99, 497 
Sheen, Y., Jeong, H., Yi, S. K., et al. 2009, AJ, 138, 1911 
Smith, B. J. & Hancock, M. 2009, AJ, 138, 130 
Smith, B. J., Struck, C, Hancock, M., et al. 2007, AJ, 133, 791 
Smith, B. J., Struck, C, Hancock, M., et al. 2008, AJ, 135, 2406 
Sulentic, J. W., Pietsch, W., & Arp, H. 1995, A&A, 298, 420 
Toomre, A. 1977, in Evolution of Galaxies and Stellar Populations, ed. 

B. M. Tinsley & R. B. Larson, 401-+ 
Weilbacher, P. M., Due, P.-A., & Fritze-v. Alvensleben, U. 2003, A&A, 397, 

545 

Weilbacher, P. M., Due, P.-A., Fritze v. Alvensleben, U., Martin, P., & 

Fricke, K. J. 2000, A&A, 358, 819 
Werk, J. K., Putman, M. E., Meurer, G. R., et al. 2008, ApJ, 678, 888 
Werk, J. K., Putman, M. E., Meurer, G. R., et al. 2010, AJ, 139, 279 
Williams, B. A., Yun, M. S., & Verdes-Montenegro, L. 2002, AJ, 123, 2417 
Wu, H., Cao, C, Hao, C.-N, et al. 2005, ApJ, 632, L79 
Wu, Y, Charmandaris, V, Hao, L., et al. 2006, ApJ, 639, 157 
Yoshida, M., Yagi, M., Okamura, S., et al. 2002, ApJ, 567, 118 
Yun, M. S., Verdes-Montenegro, L., del Olmo, A., & Perea, J. 1997, ApJ, 

475,L21 + 

Zubko, V, Dwek, E., & Arendt, R. G. 2004, ApJS, 152, 211 



